Transport through asymmetric two-lead junctions of Luttinger liquid wires 
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We calculate the conductance of a system of two spinless Luttinger liquid wires with different 
interaction strengths f/i, (?2, connected through a short junction, within the scattering state formal- 
ism. Following earlier work we formulate the problem in current algebra language, and calculate 
the scale dependent contribution to the conductance in perturbation theory keeping the leading 
universal contributions to all orders in the interaction strength. From that we derive a renormal- 
ization group (RG) equation for the conductance. The analytical solution of the RG-equation is 
discussed in dependence on gi,g2- The regions of stability of the two fixed points corresponding to 
conductance G — and G — 1, respectively, are determined. 



Introduction. Charge or spin transport in three- 
dimensional metallic systems may be described in terms 
of Landau quasi-particles. In strictly one-dimensional 
quantum wires quasi-particles are found to be unstable 
on account of the interaction between electrons. An im- 
portant part of that physics is captured by the exac tly 
solvable Tomonaga-Luttinger liquid (TLL) model.El^ As 
is well known, transport through junctions of TLL wires 
is strongly affected by the interaction in the wires, in 
some cases leading to a complete blocking of transmis- 
sion. The latter behavior can be traced back to the for- 
mation of Friedel oscillations of the charge density around 
the barrier, leading to an infinitely extended effective bar- 
rier potential in the limit of low energy. The suitable 
language to describe this situation is the renormalization 
group (RG) method, allowing to calculate the conduc- 
tance as a function of (length or energy) scale. Quite gen- 
erally, the transport behavior at low energy/temperature 
is dominated by only a few fixed points of the RG flow. 
In the neighborhood of these fixed points the conduc- 
tance is found to obey power-law behavior as a function 
of temperature T for the infinite system or as a function 
of system length L at zero temperature. The problem 
of the two-wire jun ction has been studied first in the 
pioneering worka^, using the method of bosonization. 
Later, their results have been confirmed by many other 
authors, se^ and references therein. For special values 
of the interaction exact results have been obtainecP^^. 

A purely fermionic formulation of the problem has been 
introduced in the limit of weak interaction by Yue, Glaz- 
man and Matveev^. We have extended that theory to the 
regime of strong coupling by summing up an infinite se- 
ries of terms in perturbation theory, identified as the lead- 
ing and universal contribution^^. The results obtained 
with our method for the two-lead junction as w ell as the 
three lead junction with time-reversal symmetrjl^Slill and 
with magnetic flux^^ are in agreement with exact theo- 
retical results, where available. However, our results go 
beyond what has been obtained by other authors in vari- 
ous ways. The majority of the previous works considered 



the symmetric case of equal interactio n stre ngth in both 
half wires, with exception of the work^^^^^^ discussed be- 
low. 

In this paper we generalize our previous treatment^ to 
the case of two wires with different interaction strength 
5i , 52 ■ This includes the case of a barrier at the end of 
a TLL wire, for which one of the interaction parameters 
is zero, e.g. g2 = 0. As before, we confine our considera- 
tions to spinless fermions. 

The model. To illustrate the system we are interested 
in, we first consider a tight binding Hamiltonian Jitb of 
free spin-less fermions describing two quantum wires con- 
nected at a single junction by tunneling amplitudes: 
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Here creates a fermion in wire j at site n, and ti, is 
the tunneling amplitude connecting the sites (j, n — 0) 
at the junction. The 2x2 ^-matrix characterizing the 
scattering at this junction has the structure (up to overall 
phase factors in the individual wires) 
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We choose this parametrization in terms of the trans- 
mission and reflection amplitudes i, r , since it is readily 
generalizable to the case of multi-wire junctions {n wires, 
n > 2). The above form of the 5-matrix is completely 
general (up to irrelevant phase factors) and, in fact, de- 
fines our model. Passing to the continuum limit, lineariz- 
ing the spectrum at the Fermi energy and adding forward 
scattering interaction of strength gj in wire j , we may 
write the TLL Hamiltonian % in the representation of 
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incoming and outgoing waves as 

^0 2 

H = dxY,[Hj+H}'''e{~L<x<~l)],{3) 

J-co 

= VFi^]^,jyipjAn-^F^p]^outi^^j,out, (4) 

We put vf = i from now on. The range of the interaction 
hes within the interval {l,L), where I > serves as a ul- 
traviolet cutoff and separates the domains of interaction 
and potential scattering on the junction; non-interacting 
leads correspond to large |a;| beyond L. In terms of the 
doublet of incoming fermions ^' = (V'l.m, V'2,m) the out- 
going fermion operators may be expressed with the aid 
of the S'-matrix as ^'(x) — S ■ '^{—x) . We express the 
interaction term of the Hamiltonian in terms of density 
operators pj^^n = "^^j^n'^jM = = Pj, and pj^out = 

i'j.out'^j,out = * = Pj , where pj = S+ ■ pj ■ S , a.s 

ff]"* = 2TrgipiPi + 27r52P2P2 • (6) 

The matrices are given by {pj)ap = Sa/3^aj and {fij)ap — 
S^jSjp. A convenient representation of 2x2-matrices is 
in terms of Pauli matrices aj, j = 1, 2, 3, the generators of 
SU{2) (set^). Notice that the interaction operator only 
involves us (besides the unit operator ((To)q,9 = Sap)- 
We note Tr{aj) = 0, Tria^ak) = 2Sjk, j = 0,1,2,3. 
Defining a two-component vector s = {a3,ao), we have 
Pj = ■\/l/2^^ Rj^Sfj,, where the 2x2 matrix R has the 
properties R^^ = R-'", detH = 1, and Rn — R12 = 
—R21 = R22 — l/\/2. The outgoing amplitudes will be 
expressed in terms of aj = ■ (Jj ■ S. With the aid of 
the aj the S'-matrix may be parametrized by three angu- 
lar (Euler) variables (see^), S = e*"i'^3/2g»a2<Ti/2g»a3T3/2^ 
For the case under consideration only two of these, 6, 0, 
are relevant: S = e-"^'^3/2e*('^-2e)ai/2gi0<T3/2_ -pj^g 

responding elements of the S'-matrix have been given in 
Eq. ([2|. 

Parametrization of conductance matrix. We may de- 
fine a 2 X 2 matrix of conductances Gjk relating the cur- 
rent Ij in lead j (flowing towards the junction) to the 
electrical potential Vk in lead k : Ij = X^fc^ifc^fc- I* 
follows from the conservation of charge that J^j Gjk — 
and from invariance under a shift of the zero of energy 
that Gjk — 0. Therefore, the conductance matrix 
has only one independent element G = (1 — a)/2, which 
relates the net current I — ^{Ii — I2) to the bias volt- 
age = {Vi - V2), I = GV. We note the relation 
G = i(R^GR)ii , while all other elements of R^GR 
are zero. 

In the linear response regime the conductances are 
related to the S-matrix by Gjk = Sjk — Tr(p^pk) — 



6jk — \Sjk\'^, where the label r indicates that the quan- 
tity is fully renormalized by interactions. Defining Yjk = 
Tr{pjPk) and using the above relation of Pj and s^, we 
see that the conductance components may be expressed 
in terms of F^^ = k'^riSl.s^) as Y^k = (RYR^)^.^. 
Here and in the following bold faced quantities marked 
with overbar are matrices in the transformed space (in- 
dices p,iy )■ It follows from the analysis of the trans- 
formed matrix R-^GR given above that the matrix Y 
has block structure, the nonzero elements being given 
by the conductance parameter introduced above and by 
unity, Yii = a, Y22 = 1, Y12 = Y21 = 0. From the 
above relations we see that the parameter a may be ex- 
pressed by the angle 9 in the above parametrization of 
the S-matrix as a = — cos 20. We find therefore that a is 
confined within the region a € [—1, 1]. 

Perturbation theory. We now calculate the conduc- 
tance in perturbation theory in the interaction. In 
first order we have to evaluate the diagrams depicted 
in Fig. [1] Here solid lines denote Green's functions in 
position-energy representation (a;„ are fermionic Matsub- 
ara frequencies), Q{x]ujn) = —isign{ujn)0{xujn)e~''^'^". 
The double wavy lines denote the interaction opera- 
tor, which in the lowest order is given by the matrix 
g^^ = (R^gR)^^, with gjk = Sjtgj, in the transformed 
or initial representations, respectively. At the ends of 
the interaction lines operators (73 and (73 are attached 
depending on whether a; < or a; > 0. As a result, one 
finds in lowest order in the interaction 



= -^Sf^iS^iia^ - l){gi+g2)A. 

Rerejhe W^, = (R^WR)^^ and W,k = 

[pj,Pk] , are 2x2 matrices for each pair of pv (or jk) and 
the trace operation Tr is defined with respect to that 
matrix space, whereas the g^^ = (R^gR)^^, with gjk = 
Sjkgj are scalars. Notice that only one matrix element, 

Wii, is nonzero. 

We may extend this analysis into the strong coupling 
regime by summing up infinite classes of contributions in 
perturbation theory. It has been found in^ that the dia- 
grams shown in Figs. [T] and [2] provide the leading terms in 
the neighborhood of the fixed points. Their contribution 
is universal in contrast to other higher order terms (see 
below). These may be interpreted as a renormalization of 
the bare interaction, 2-!TgkiS(x — y) — ^ Lkiix,y;u!n). The 
contribution to conductance stemming from the first two 
diagrams (a) and (b) in Fig. [T] in this ladder approxima- 
tion is given by 
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FIG. 1: Feynman diagrams depicting the corrections to conductance. Two first diagrams, (a) and (b), correspond to Eq. ([8|, 
the third diagram ("vertex correction") vanishes in the static limit, external Q — > 0. 



G("+^) = --T^ V / dxidx2 dyLn{xi,X2;uj)g{y + x-e + n) 

X {g{xi - x; e)gi-xi - X2;e~ uj)g{x2 - y; e) + cos 46* tj(-a;2 - x; e)g{xi +X2;e + uj)g{~xi -y;e)) 



(8) 



Here only the (1, 1) element of L enters, corresponding to the fact that operators CT3 and CT3 are attached to the ends 
of the renormalized interaction line. The factor cos 40 appears as Ti[a^a'iaza^) /2. In the limit of zero temperature 
we may convert the summation over Matsubara frequencies to an integration along the imaginary axis. Another 
contribution to G is obtained by reverting the arrows on the fermion lines in Fig. [T^, Fig. [TJd, and doubles the above 
result. Performing the integrations on e and on y from L to 00 and taking the limit f2 — >■ we find 



G^^' = (1 — a^) / dxidx2— Lii{xi,X2',u})9{uj)e 

ZTT 



It is useful to first calculate L in the initial representation, where the interaction matrix g is diagonal, but the matrix 
Y is nondiagonal. Then L is found to satisfy the following linear integral equation for > (Ref. [5], note that the 
definitions of Y and for the strength of interaction there are different): 



L(a;,y;w) 
L2(a;,2/;a;) 



2n^5{x - y) ( J ) - 2^ / dz 



gYn(x + z,w), ^Ii{x- z,ijj)\ fL{z,y:uj) 



gn(z - x,uj), 







L2(^,2/;w) 



with the fermionic loop Il{x,u!n) = (Stt) ^{S{x) — \uJn\0{xuJn)e ^""). Expressing 'L2{x,y;uj) via L(a;,y;a;) we have 



L{x, y; io) = 2ngS{x - y) + ojg j dz {Y - ig)e-"(-+^) - ige-^l^-^l 



with gj ~ gj/d'j and dj = y 1 ^ 9j- order to solve this integral equation we first calculate a partial summation, 
defined by C{x,y;uj) = limy_j.o L(a;, y; w). Then L(x,y;a;) will be the solution of 



.{x, y; uj) = C{x, y;uj) + -^ / dzidz2C{x, zi; w)Ye-"(^i+^^)L(z2, y; cj) 



(9) 



and C{x,y;uj) satisfies the integral equation 



C(a;,y;a;) = 27rg5(a; - y) - iwgg / dz 



C{z,y;uj) . 



(10) 



These integral equations are shown diagrammatically in Fig. [2] 
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We now define C{x,y;uj) = 2TTg5{x — y) — Ci{x, y;uj), so that the inhomogeneity in the integral equation for 
Ci{x, y;uj) is differentiable 



Ci{x,y;oj) = ttujI 



The integral equation for Ci(x,2/;a;) may be converted into a second order differential equation 




FIG. 2: Feynman diagrams depicting the integral equations 
for the renormalized interaction, Eqs.([9]l and (10 1. 



dx'^ 



ggj 



Ci(a;,y;a;) = -27ra; g gj(x- y) 



Since the matrix g is diagonal, Ci is diagonal and 
we have two uncoupled differential equations for the 
components Ci^-, j = 1,2. Taking into account the 
boundedness of Cij-, the general solution is given by 
Cij(a;,2/;cj) = ^^jg,g]{Aj[ij)e-^^'' + e-^^^-v\) , where 
ijjj = uj/dj. It follows from the boundary conditions at 
a; = that Aj{y) = e~^^^ . The quantity C is thus a 
diagonal matrix given by 



C(x,2/;w) = 27rd '^^5{x-y) 



TTCjd ^g^ 



We now return to the integral equation for L(a;,2/;a;). 
Since the kernel is separable it is easily solved. 
We define the auxiliary matrix functions U(?/;a;) = 
/(ia;e-'^^L(x,?/;w), V(y;w) = J dxe-'^'=C{x,y;uj) = 
y^'^ as well as the matrix = 
-"(="+f)C(a;,y;w) =g(l + d)-i . Multiplying 
the above integral equation ^ by e^"^^ and integrating 
over x we find U = V + Q-^YU = (1 - Q-^ Y)-iV and 
we finally get 



27rd-ige-" 
5¥ / '^^dye 



2tt 



Yix;Lu)Yil-Q-'Y)-'Viy;u) 



This result is now substituted into the expression for 
the conductance. Performing the integration over y 
first we get / dye-'^yL{x,y;uj) = Y{x;uj)[l + Y(l - 
Q-iY)-iQ-i] = V(a;;a;)[l - YQ-i]-i . Next inte- 
grating over positive w yields / (da;/27r)e^"^V(x;a;)[l — 
YQ-y = (l/x)g(l + d)-i[l-YQ-Vi = (l/x)[Q- 
Y] ^. Finally, the integration over x produces the scale 
dependent logarithm A = ln{L/l) . The conductance 
in the ladder approximation is found after transforming 
to the rotated basis (quantities with overbar) and taking 
the (1, l)-element of the matrix. 



G(^)^i(l-a)-(l-a2) [(Q-Y)-i]^^A 

Renormalization group equations. The renormalization 
of the conductances by the interaction is determined from 
the scale dependent contributions in perturbation theory. 
Differentiating these results with respect to A (and then 
putting A = 0) we find the RG equation for the quantity 
a = 1 — 2G in the ladder approximation 



da 
dK 



2(l-a2)[(Q-Y)-i]^^ (11) 



Here Y — diag{a, 1) and Q = R-^QR , such that Qn — 
Q22 = Q+ = ill + 92)/2 and Q^^^ = Q21 = Q- = (Qi - 
92)72 where q, = (1 + ^1 - g])/9, - (1 + K,)/{1 - 

Kj), and Kj ~ ^(1 — gj)/{l + gj) is the usual Luttinger 
liquid parameter for wire j. We define 



7-Q+-0-/(Q+-i) 



(12) 



(note that I7I > 1 for any Ki_2 > ), then the RG- 
equation takes the explicit form 



da ^o? 
dK a 



- 1 



■ 7 



(13) 



RG flow and conductance. The fixed points of the 
above RG equation are labelled , at a = 1, G = 0, 
(complete separation of the wires) and A at a = — 1, 
G = 1 (ideal conductance through the junction). In order 
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FIG. 3; Stability regions, labeled by the corresponding fixed 
point are shown in K1-K2 plane. The hyperbola separating 
the regions is shown by solid line, the non-interacting values 
i^i,2 = 1 are shown by dashed lines. 



to discuss the stability and to calculate the conductance 
we rewrite the RG-equation in terms of the conductance 



dG_ Gjl -G) _ 



(14) 



Stability of a fixed point requires that the derivative of 
/3(G) at the fixed point is negative. At fixed point N this 
translates into 7 — 1 > 0, or more explicitly, qi + (?2 > 2, 
and in terms of the Luttinger parameters K^^ +K2^ > 2 
i.e. either both interactions should be repulsive, or one 
of them is attractive, but weak relative to the repulsive 
one. At fixed point A the condition is 7 + 1 < 0, or 
Q1+Q2 < 2, and therefore K^^ +K2^ < 2 , meaning that 
the interaction is predominantly attractive, but a weaker 
repulsive interaction in one of the wires is possible. The 
line separating the stability regions in the ifi-ii'2-plane 
is the hyperbola, {Ki — ^){K2 — ^) = |, passing through 
the (no interaction) point Ki — K2 — 1- This hyperbola 
and corresponding stability regions are shown in Fig. |3] 
The RG-equation may be integrated to give 



G^~^[1~GY+'' = c{L/l)^ (15) 

In the vicinity of fixed points N and A , respectively, we 
then find the power laws 



G = cn{IIlY'^^~^\ at TV (16) 
G = l-CA{L/lf'^^+^\ at A (17) 

Exphcitly we have 4/(7 - 1) = (i^f^ + K^^ - 2), and 
4/(7 + 1) = 2{Ki +K2- 2KiK2)/iKi + K2). At Ki = 



K2 = K we recover the well- knowrP^ exponents 2{K ^ — 
1), 2(1-/^).^ 

Of special interest is the case 171 — g, 92 — 0, when 
we have 7 = 2g - 1 = (1 + 3K)/{1 - K) and the 

power laws are given by G = cn{L/1)^''^ at N and 
1 - G = CA(i/Z)-2(^-i)/('^+i) at A. We may compare 
these findings with results obtained in Ref. |13j for a non- 
symmetric position of an impurity in a TLL wire. When 
taking the impurity position sufficiently close to the in- 
terface between the TLL wire and a non-interacting lead, 
scaling exponents identical to the above ones were found 
irfi^, as can be read off from Eqs. (10) and (13) there. 

Thus we recover the correct exponents in the scaling 
law for the conductance. We obtained them by summa- 
tion of the ladder sequence for the renormalized interac- 
tion in the presence of the junction. As was shown irP, 
the first contributions to conductance beyond the ladder 
series appear in the third order of interaction and are 
of three-loop type (the two-loop RG contributions are 
absent). These latter terms do not contribute to above 
scaling exponents, but define the relation between the 
prefactors cm and ca- The Eq. (151 implies that both cn 



and Ca depend on the initial conditions, encoded in the 
amplitude c. It is easily seen, that the ladder approxi- 
mation corresponds to the relation 



Vci+^ = 1 



It was showrP that the three-loop corrections change this 
ratio to values, different from unity and depending on 
the strength of interaction and regularization scheme, i.e. 
non-universal. 



Our expression ( 12 ) shows that the boundary expo- 
nents depend only on the sum K^^ -\-K2^ ■ This is in pre- 
cise agreement with Eq. (13) in ReLlSl^ where the scaling 
exponent of the point contact between two chiral (Hall 
edge) states was derived. On the basis of this observa- 
tion it was suggested there, that the combination of two 
chiral states, Ki = 1/3 and K2 = 1, is equivalent to 
the well-known problem of impurity in TLL wire with 
Ki = K2 = 1/2, which can be fully solved.'^ We note 
here that the coincidence of the scaling exponents at two 
fixed points N and A might not necessarily mean the 
coincidence of the full scaling form for the conductance. 
The non-universal three-loop terms (which are not dis- 
cussed in the standard bosonization approach) may be 
different in cases Ki = 1/3, K2 — 1 and Ki = K2 = 1/2. 

Conclusion. In this paper we employed a fermionic de- 
scription of a general two-wire junction of two TLL-wires 
to derive the renormalization group equation for the con- 
ductance, using the approach developed by us earlieJ^. 
We used an infinite summation of perturbation theory 
in the form of a ladder approximation, allowing for an 
analytical solution for arbitrary junction parameters and 
interactions in the wires. As demonstrated earlier the 
approximation is asymptotically exact in the vicinity of 
the fixed points. As in the well-studied case of a sym- 
metric junction there exist two fixed points of the RG 
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flow. Fixed point N corresponds to a complete separa- 
tion of the wires, i.e. the conductance vanishes. It is 
stable in a region of the A'i-/f2— coupling constant plane 
which is predominantly repulsive, meaning that a weakly 
attractive component, say K2 > 1 is permitted. At fixed 
point A the conductance assumes its maximum value. 
It is stable in the complementary part of the coupling 
constant plane. The two stability regions in the Ki- 
K2— plane are separated by a hyperbolic boundary curve. 
The representation chosen is maximally general and may 
therefore be easily extended to junctions connecting more 



than two wires. The 1-2-symmetric three wire junction 
( "Y-junction" ) has been considered by us in Refs.tl™!^'^. 
Work on the four- wire junction is in progress. 
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